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Angular distributions of elastic 140 Ce(a:,a) 140 Ce scattering are analyzed in the framework of 
the optical model from low energies around the Coulomb barrier up to about 40MeV. From the 
local fits the total reaction cross section <7 rcac is extracted. This procedure requires experimental 
scattering data with small uncertainties. The results for cr reac are compared to predictions of different 
systematic global potentials. It is found that the total reaction cross section cr reac is well predicted 
from all global potentials under study although the reproduction of the angular distributions is 
not perfect in all cases. In addition, the lower energy limit for the extraction of <T roac from elastic 
scattering angular distributions is analyzed. Finally, the potentials under study are used to calculate 
the 143 Nd(n,a) 140 Ce cross section, and the predictions are compared to experimental data. 

PACS numbers: 24.10.Ht,24.60.Dr,25.55.-e,25.55.Ci 



I. INTRODUCTION 

The total reaction cross section <7 re ac is a basic quan- 
tity for the prediction of reaction cross sections in the 
statistical model. It has been found that especially at 
low energies huge discrepancies are found between predic- 
tions for a-induced reaction cross sections using different 
global a-nucleus optical potentials. This holds in partic- 
ular for (a/y) capture reactions for targets with masses 
above A m 100 [lHg], but also the energy dependence of 
recent (a,n) data for 141 Pr was difficult to fit [jjl [To|. 

The total reaction cross section cr reac is related to the 
complex scattering matrix Sl — Vl exp (2££r,) by the 
well-known relation 

a rC ac = ^^ = ^^(2i + l)(l-r / !) (1) 

L L 

where k = y/2^,E c , m Jh is the wave number, E c m . is the 
energy in the center-of-mass (cm.) system, and r\i, and 
Sl are the real reflexion coefficients and scattering phase 
shifts. &L is the contribution of the L-th partial wave 
to the total reaction cross section er reac . The relation in 
Eq. ([1]) has recently been confirmed experimentally at 
low energies fill ]. 

The present study determines a reac from the angular 
distributions at E lah = 14.7 MeV and 19.0, 24.0, 
32.0, and 37.7 MeV [Hj]. The latter data are unfortu- 
nately not published. The data have been measured at 
the XTU Tandem of the INFN Laboratori Nazionali di 
Lcgnaro using a cerium target enriched to about 96 % 
in the semi-magic N = 82 nucleus 140 Ce. Access to the 
numerical 140 Ce(a,a) 140 Ce scattering data of [l3| will be 
provided via the EXFOR database [14| • In addition, elas- 
tic scattering data for natural cerium ( 140 Ce: 88.48%) 
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are available in literature. These data have also been 
analyzed, and the results fit into the observed systemat- 
ics of total reaction cross sections. However, the angular 
distributions are not shown here because they may be 
affected by the other cerium isotopes. 

The obtained results for £r rC ac are compared to predic- 
tions from recently published a-nucleus potentials (l6r - 
[H| and to the widely used potential by McFadden and 
Satchler The example of 140 Ce is well-suited for a 
study of uncertainties for £r rC ac at low energies because 
experimental scattering data are available down to rela- 
tively low energies. Further information on the a-nucleus 
potential can be derived from reaction data at very low 
energies. Often the M reaction has been used for this 
purpose (see e.g. H flo|). but in the present case there 
are no experimental data for the 140 Ce(a,n) 143 Nd reac- 
tion. Instead, the reverse 143 Nd(n,a) 140 Ce reaction will 
be investigated here. 



II. OPTICAL MODEL ANALYSIS 

The elastic scattering angular distributions are ana- 
lyzed in the usual way within the framework of the op- 
tical model and a complex a-nucleus optical potential. 
Details of the fitting procedure for this local fit are simi- 
lar to Refs. pilllll and are only briefly summarized here. 

The total potential U(r) is composed of the real 
Coulomb potential Vc (V) and the complex nuclear poten- 
tial V(r) + iW(r). The Coulomb potential is calculated 
from a homogenously charged sphere where the Coulomb 
radius is taken from the root-mean-square (r.m.s.) radius 
of the real folding potential (see below). Properties of the 
real and imaginary part of the nuclear potential are pre- 
sented in the next paragraphs. 
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A. Real part of the nuclear potential 

The real part of the nuclear potential is determined by 
a double-folding procedure where the folding potential 
Vp (r) is modified by a strength parameter A and a width 
parameter w: 



V(r) =XV F (r/w) 



(2) 



The strength parameter A and the width parameter w 
will be adjusted to the experimental 140 Ce(a,a) 140 Ce an- 
gular distributions. Obviously, the width parameter w 
should remain close to unity; otherwise, the folding po- 
tential would be questionable. The strength parameter A 
is typically around 1.1 to 1.4 leading to volume integrals 
per interacting nucleon pair of Jr « 310 — 350MeVfm 3 
[20(. (As usual, the negative signs of Jr and Jj are ne- 
glected in the following discussion.) 

The calculation of the folding potential requires the 
density of the 140 Ce target nucleus which is usually de- 
rived from electron scattering. Although no data are 
available in the widely used compilation 21], two pa- 
pers [13, HI] use a 140 Ce density and refer to a "com- 
plementary, high-precision, clastic electron scattering ex- 
periment performed at Saclay" with the reference "D. 
Goutte et al, to be published". The density parameters 
in Fourier-Bessel parametrization can be found in an un- 
derlying Ph.D. thesis by B. L. Miller [H, Table A.l, and 
are repeated here for easier access in future work (see 
Table |T|. The resulting folding potential is very close to 
the average of the potentials for the neighboring N = 82 
nuclei 139 La and 142 Nd; thus, the reliability of the unpub- 
lished electron density of 140 Ce from [24| is confirmed. 



TABLE I: Experimental charge density distribution of 140 Ce 
as listed in Table A.l of [U, parametrized by a Fourier-Bessel 
series: p(r) = Gn x jo(nirr/RFB) with the Fourier-Bessel 
cutoff radius Rfb = 10 fm. 



Gn 



1 0.759956 2 0.578531 

3 -0.589376 4 -0.184926 

5 0.337939 6 0.006674 

7 -0.132665 8 0.002477 

9 0.002589 10 -0.013456 



B. Imaginary part of the nuclear potential 

The imaginary part of the optical potential is 
parametrized by the usual Woods-Saxon potentials of 
volume and surface type. Except at the highest energy 
under study where an additional volume Woods-Saxon 
potential is required, only a surface Woods-Saxon poten- 
tial was used: 



W(r) = W v f(x v ) + W s 



df{x s ) 
dxs 



Wi are the depth parameters of the volume and sur- 
face imaginary potential, and the Woods-Saxon function 
f(xi) is given by 



f(xi) = 1 + cxp (Xi) 



(4) 



with Xi = (r — RiAj, )/a% and i = V, S for the volume 
and surface part. Note that Wy < and Ws > in the 
chosen conventions Q and (H} for an absorptive negative 
W(r) < 0. The maximum depth of the surface imaginary 
potential is given by — Ws/4 at r 



i?s4 /3 - 



Results 



(3) 



The parameters of the nuclear potential are adjusted 
to the experimental angular distributions using a stan- 
dard x 2 minimizing procedure. The resulting parameters 
are listed in Table HH and the fits are compared to the 
experimental angular distributions in Fig. [1] Excellent 
agreement between the local fits and the experimental 
data is obtained at all energies under study. For the real 
part the parameters of the fits show very minor varia- 
tions with energy and are in their expected ranges [20| . 
The parameters of the real and imaginary parts of the 
potential will be discussed in detail below. 

The width parameter w is about 1.5% larger than 
unity and practically constant except at the highest en- 
ergy for 140 Ce. The strength parameter A varies by about 
10 % between 1.15 and 1.26; this leads to a similar varia- 
tion of the real volume integral Jr which shows a maxi- 
mum of 341 MeVfm 3 at 24MeV and slightly smaller val- 
ues at higher and lower energies (similar to the finding 
in [HI)- The variation of the potential parameters of the 
real part is relatively small, and thus the real part of 
the potential and its energy dependence should not be 
the reason for major uncertainties in the prediction of 
a- induced reaction cross sections at low energies. 

Contrary to the real part, the imaginary part shows a 
significant energy dependence. The most obvious signa- 
ture is the change from a surface Woods-Saxon poten- 
tial at very low energies to a dominating volume Woods- 
Saxon potential at the higher energies under study. It 
is sufficient to use a pure surface Woods-Saxon potential 
at the lower energies below 25MeV, and reduced x 2 /F 
values around unity are obtained. The variation of the 
geometry of the surface imaginary part at the lower en- 
ergies is small: the radius parameter Rg varies by less 
than 5 % around its average value, and the diffuseness as 
shows a somewhat larger spread of about 25 % around its 
average. However, there is no systematic energy depen- 
dence of Rs and as which may lead to uncertainties in 
the extrapolation of the potential down to very low ener- 
gies (see also Sect. IIII B]) . The depth parameter Ws and 
the imaginary volume integral Jj increase with energy 
for the lowest 3 energies. 

At higher energies above 30MeV, the best x 2 /F is ob- 
tained from a combination of a volume and a surface 
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TABLE II: Parameters of the optical potential and the total reaction cross section a rcac derived from 140 Ce(a,a) 140 Ce angular 
distributions. 



Slab 


E c .m. 


A 


w 


Jr 




W v 


Rv 


By 


W s 


Rs 


as 


Ji 


^/,rras 


X 2 /F 


a 

Oreac 


Ref. 


(MeV) (MeV) 


(-) 


(") 


(MeVfm 3 ) 


(fm) 


(MeV) 


(fm) 


(fm) 


(MeV) 


(fm) 


(fm) 


(MeVfm 3 ) 


(fm) 


H 


(mb) 


Exp. 


14.69 


14.28 


1.146 


1.021 


315.8 


5.641 








28.0 


1.444 


0.507 


18.2 


7.774 


0.1 


28±4 


[12] 


19.0 


18.47 


1.218 


1.014 


328.1 


5.599 








101.8 


1.517 


0.303 


43.1 


7.970 


0.7 


399±30 


[13] 


23.97 


23.30 


1.253 


1.017 


340.7 


5.616 








114.0 


1.445 


0.434 


63.2 


7.709 


0.8 


999±30 


[13] 


32.0 


31.11 


1.276 


1.006 


336.0 


5.557 


-21.1 


1.311 


0.211 


13.4 


1.574 


0.416 


58.6 


5.860 


0.8 


1453±44 


[13] 


37.7 


36.65 


1.265 


0.999 


326.3 


5.518 


-19.1 


1.374 


0.269 


8.9 


1.580 


0.729 


62.6 


6.218 


1.0 


1824±55 


[13] 


45.0 6 


43.75 


1.195 


1.020 


328.3 


5.635 


-28.2 


1.028 


0.809 


28.0 


1.516 


0.504 


59.3 


6.291 


0.4 


1940±97 c 


[15] 



a from the local potential fit using Eq. (TTJ; uncertainties estimated 
from the model-independent phase shift analysis 
''natural cerium ( 140 Ce: 88.48%) 

c increased 5 % uncertainty because of contributions from other 
cerium isotopes 



Woods-Saxon potential. At all energies above 30 MeV 
the volume term is dominating. An interesting ambigu- 
ity in the parametrization of the imaginary potential is 
found for the combination of volume and surface Woods- 
Saxon poentials. Very similar potentials can be obtained 
using a positive or a negative surface term. This is il- 
lustrated in Fig. [2] for the analysis at E = 31.11 MeV. 
In the first fit (parameters given in Table [TTJ) the usual 
behavior is found: the weak negative surface component 
modifies the dominating volume component in the sur- 
face region around 8 fm. However, in the second fit (pa- 
rameters: A = 1.269, w = 1.004, W v = -23.1 MeV, 
R v = 1.579 fm, a v = 0.525 fm, W s = -64.6 MeV, 
R s = 1.394 fm, a s = 0.571 fm, cr reac = 1464 mb) a rela- 
tively strong and positive surface contribution is found. 
But the total imaginary part is very similar in both fits, 
in particular in the most relevant surface region, and 
thus also the resulting jc/F and derived cx reac are practi- 
cally identical. Such an ambiguity is found at all energies 
above 30 MeV. In principle, the fit with the lowest x 2 /F 
should be used for the extraction of £r rC ac from the experi- 
mental angular distribution. However, because the found 
X 2 / F minima arc very similar and do not differ by more 
than 10 %, and because the derived <7 rea c are stable within 
1 %, the fits with the usual negative surface component 
are listed in Table [TTJ This choice should avoid further 
complications e.g. when using these parameters for the 
construction of a global potential. For completeness it is 
pointed out that fits with a pure volume Woods-Saxon 
imaginary potential have a worse x 2 1 F, and the derived 
froac vary within several per cent. 

It is obvious that the energy dependence of J/ and the 
shape of the imaginary part of the potential are dominat- 
ing sources for uncertainties of cross section predictions 
at very low energies. As usual, the imaginary potential 
W(r) is dominated by the surface contribution at low en- 
ergies, and its volume integral Jj increases with energy. 
At higher energies above 30 MeV the volume contribu- 
tion becomes dominant. The data at 23.30 MeV could 



also be fitted using a combination of volume (Wv = 
-13.8MeV, R v = 1.30fm, a v = 0.81 fm) and surface 
(W s = 13.6 MeV, R s = 1.58 fm, a s = 0.34 fm) poten- 
tial, i.e. with a further decreasing volume contribution. 
Although the number of ajustable parameters is larger, 
this fit did not lead to an improved description of the an- 
gular distribution (no significant improvement of \ 2 I F). 
Therefore, from this energy down to lowest energies only 
a surface imaginary potential was used. This choice is 
confirmed by the fact that an excellent reproduction of 
all angular distributions below 25 MeV with x 2 /F < 1 
has been achieved using a pure surface Woods-Saxon po- 
tential in the imaginary part. Because the volume in- 
tegrals J j increase with energy at these low energies, a 
similar behavior is found for the strength of the surface 
imaginary potential Ws- 



The influence of uncertainties in the real part and the 
imaginary part on the calculation of the total reaction 
cross section <7 roac will be investigated in further detail 
later (see Sect. IIIip . 



The elastic scattering angular distributions are also 
compared to predictions from several global a-nuclcus 
optical potentials. In Fig. [JJ we show the very simple 4- 
mrameter potential by McFadden and Satchler (MCF) 
IS | , the many-parameter potential by Avrigeanu (AVR) 
16(, and the recent first version of the ATOMKI potential 
(ATOMKI-V1) [Hj]. For illustration, I also show the real 
and imaginary nuclear potentials for the lowest energy 
of 14.28 MeV and an intermediate energy (31.11 MeV) in 
Fig. El 



4 



1.1 
1.0 

as 
b 

b 0.9 
0.8 



10 J 



a) 



experimental data: 
B. D. Watson et at, 
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FIG. 1: (Color online) Rutherford normalized elastic scat- 
tering cross sections of the 140 Ce(a,a) 140 Ce reaction versus 
the angle in center-of-mass frame. The lines are calculated 
from a local potential fit which is adjusted to the scatter- 
ing data (full black line) and from different global a-nucleus 
potentials [H, El El (MCF: blue dotted; AVR: red dash- 
dotted; ATOMKI-V1: green dashed). The experimental data 
are taken from [l^. Il3j|. The given energies are E c , m . in the 
center-of-mass system. 



III. DETERMINATION OF THE TOTAL 
REACTION CROSS SECTION a rcac AND THE 
REDUCED CROSS SECTION a Icd 

A. Extraction of d rcac and ovcd from the 
experimental angular distributions 

The total reaction cross sections cr rcac are derived from 
these local fits using Eq. ([1]) . For an estimation of the un- 
certainty of (Treac an additional model-independent phase 
shift analysis (PSA) has been performed using the tech- 
nique of [261 ] . The total reaction cross sections <7 rea c from 
the PSA are close to the results of the local potential fit. 
A typical uncertainty for cr roac of about 3 % is estimated 




r(fm) 

FIG. 2: (Color online) Ambiguity of the imaginary poten- 
tial, illustrated for the analysis of the angular distribution at 
E — 31.11 MeV. The imaginary potentials of two different fits 
(shown in blue and red full lines) are very similar. However, 
the decomposition into the volume (thin dashed) and surface 
(thin dotted) components shows that in one case there is a 
weak negative surface term (red) whereas in the other case a 
strong positive surface term (blue) is found. 



(see also [27|) except at very low energies significantly 
below the Coulomb barrier. 

It is interesting to note that there is only a relatively 
weak dependence of the derived total reaction cross sec- 
tion tJreac on the absolute normalization of the elastic 
scattering angular distribution. At relatively low energies 
the absolute normalization is well-defined from Ruther- 
ford scattering at forward angles. This is the case for 
all energies under study in this work. But even at en- 
ergies far above the Coulomb barrier it turns out that 
the potentials are relatively well-defined by the diffrac- 
tion pattern in the angular distribution. Consequently 
the extracted <7 roac vary typically by less than 10 % even 
for strong variations of the absolute normalization of up 
to 30%. 

In general, the extraction of the total reaction cross 
section cr reac from clastic scattering angular distributions 
requires theoretical considerations and is thus somewhat 
model-dependent. But because of the small sensitivity 
of fj rca c to the chosen model, the total reaction cross 
section <7 reac can be considered as a quasi-experimental 
quantity. This holds in particular for the present case 
where the elastic scattering angular distributions cover a 
wide angular range and the scattering energies are not 
extremely low. These quasi-experimental total reaction 
cross sections are compared to predictions of global a- 
nucleus potentials in Table IIIII The very simple MCF 
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FIG. 3: (Color online) Real part (a,c) and imaginary part 
(b,d) of the underlying potentials for the angular distribu- 
tions at 14.28 MeV (a,b) and 31.11 MeV (c,d) from the lo- 
cal fit (black full line) and the global potentials [l|| [HI . [IsJ 
(MCF: blue dotted; AVR: red dash-dotted; ATOMKI-V1: 
green dashed; colors and linestyles identical to Fig. [T]). A 
decomposition of the best-fit potential at 31.11 MeV into its 
volume and surface components is also given (thin dashed and 
dotted black lines). 



potential [19j, the many-parameter AVR potential |16| . 
and the few-parameter ATOMKI-V1 potential [l8[ pro- 
vide almost identical cr r cac in good agreement with the 
experimental results although the reproduction of the 
experimental angular distributions is not perfect in some 
cases (see Fig. [I). Obviously, the global potentials (which 
are not fitted to the data) cannot reproduce the experi- 
mental angular distributions with the same quality as the 
local fit. In particular, the ATOMKI-V1 potential has 



been adjusted to low-energy scattering data only. This 
leads to a relatively poor description of the experimental 
angular distributions at higher energies because of the 
missing volume term in the imaginary part. However, 
even deviations in the angular distribution of up to one 
order of magnitude at very backward angles do not pre- 
vent the prediction of cr rea c with minor uncertainties. 



TABLE III: Experimental total reaction cross sections <T roac 
(in mb) derived from 140 Ce(a,a) 140 Ce elastic scattering an- 
gular distributions (from Table [TTjl . compared to predictions 
from the global a- nucleus potentials MCF [ljj, AVR [l(|, and 
ATOMKI-V1 QJ]. For a discussion of the 143 Nd(n,a) 140 Ce 
data see Sect. IIVI 



Bc.m. (MeV) 


exp. 


MCF 


AVR 


ATOMKI-V1 


10.0 a 


(1.9 ±0.4 jib) 6 


9.3 /lib 0.7 Mb 


4.5 ^b 


14.28 


28±4 


26 


23 


24 


18.47 


399±30 


464 


477 


460 


23.30 


999±30 


987 


1025 


989 


31.11 


1453T44 


1472 


1508 


1494 


36.65 


1824T55 


1676 


1712 


1712 


43.75 c 


1940T97 


1852 


1883 


1903 



"corresponds to the low-energy 143 Nd(n,a) 140 Ce data 
''derived from the low-energy 143 Nd(n,a) 140 Ce data 
c natural cerium ( 140 Ce: 88.48%) 



The similarity of the total reaction cross sections <r reac 
from the different global potentials can be understood 
from a semi-classical interpretation of the reflexion co- 
efficients tjl in the partial wave analysis. The following 
discussion will focus on the shown energies of 31.11 MeV 
above the Coulomb barrier and 14.28 MeV below the 
Coulomb barrier (see Fig.|3]for the potentials and Fig. 0] 
for the reflexion coefficients tjl and the contributions <jl 
of the L-th partial wave to cr reac ). The discussion extends 
a similar study in [Io| . 

Small angular momenta L correspond to small impact 
parameters or central collisions. At energies above the 
Coulomb barrier the incoming a particles are able to 
reach the interior, and the corresponding partial waves 
are almost fully absorbed: j]l ps for L < 10 and thus 
cl ~ 0L m ax = (2L + \)tt / k 2 . The linear increase of is 
nicely seen in Fig. U (indicated by the brown line) . Con- 
sequently, only very few a particles are scattered into the 
backward angular region where the cross section is found 
to be much smaller than the Rutherford cross section. 

Large angular momenta L correspond to large impact 
parameters or peripheral collisions and scattering to for- 
ward angles. These partial waves do not reach the nuclear 
interior and are thus practically not absorbed: tjl ~ 1 
and <jl ~ is found for L > 20. 

Potentials have to fulfill two simple criteria to provide 
Crcac above the Coulomb barrier with sufficient accuracy. 
First, the potential must have the correct short range; 
then t]l ~ 1 for L > 20 is automatically obtained. Sec- 
ond, the imaginary part must be sufficiently strong to 
guarantee r\h ~ for L < 10. These two simple cri- 
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FIG. 4: (Color online) Reflexion coefficients r\h (a,c) and 
contributions <tl of the L-th partial wave to the total reac- 
tion cross section a rcac (b,d) for the angular distributions at 
14.28 MeV (a,b) and 31.11 MeV (c,d). The data are obtained 
from the local fit (full black squares) and from the global po- 
tentials under study HE EE El (MCF: blue circles; AVR: red 
triangles; ATOMKI-V1: green diamonds; colors are identical 
to Fig. [TJ. The data points are connected by thin colored 
lines to guide the eye. The maximum values 

f^L.max are in- 
dicated as thick brown line (31.11 MeV, d) and as data point 
for Co, max (14.28 MeV, b). Further discussion see text. 



teria are necessary but not yet sufficient for an excellent 
a-nucleus potential which should simultaneously describe 
elastic angular distributions with reasonable x 2 / F. Both 
simple criteria are fulfilled by all global potentials in this 
study. The origin for differences in cr roac from different 
realistic potentials can only stem from a few r\h orur, for 
intermediate angular momenta 10 < L < 20. But even in 
this intermediate L range the ol have to decrease from 



Q L,max down to ctl ~ 0. Thus, it is not at all surpris- 
ing that different global potentials provide very similar 
CT r0 ac at energies above the Coulomb barrier. As a test I 
have increased the imaginary part of the local potential 
at 31.11 MeV by a factor of 2 and a factor of 10 (!). As 
expected, this leads to only very moderate enhancements 
of the total reaction cross section (T reac by about 1 % and 
8 % although the corresponding clastic angular distribu- 
tions change by more than one order of magnitude in the 
backward angular region. 

The situation is much different at lower energies below 
the Coulomb barrier. Again, partial waves with large 
angular momenta L are not much affected by any short- 
range potential, and tjl ~ 1 is found. However, even a 
particles with small impact parameters will mostly not 
be able to tunnel through the Coulomb barrier. Thus, 
the elastic cross section at backward angles rises and ap- 
proaches the Rutherford cross section, and the rj^ re- 
main close to unity even for small L. The total cross 
section er rca c is now composed of the contributions from 
a few partial waves with low L < 8 in the chosen exam- 
ple of 14.28 MeV. It is surprising that even under these 
conditions the various potentials under study provide 
very similar er rca c (see Table [TTTJ) , and indeed significant 
variations have been found in the analysis of the recent 
141 Pr(ev,n) 144 Pm reaction data at slightly lower energies 



B. Sensitivity of a lea . c to details of the optical 
potential 

Let me now artificially decompose the absorption into 
an interior contribution (r < R a + RnoQ e ) and an exte- 
rior contribution (r 3> R a + Rnoc e ). It is the aim of this 
analysis to find out the radial range of the potential which 
mainly defines <x re ac- A semi-classical decomposition as 
e.g. presented in (28j may be misleading here because 
the semi-classical approximation deviates from a fully 
quantum-mechanical calculation by several per cent [28| 
which is not acceptable at very low energies. For com- 
pleteness it has also to be pointed out that "exterior" in 
the above sense means radii above sa (10 — 12) fm where 
the nuclear real and imaginary potentials have practi- 
cally dropped to values close to zero (see Fig. E]). Conse- 
quently, "interior" corresponds to the radial range where 
the nuclear potentials deviate from zero, and "interior 
absorption" is dominated by the surface absorption by 
the surface imaginary potential. 

The interior contribution requires that the incoming a- 
particlc tunnels through the Coulomb barrier and reaches 
an area where it is absorbed by the imaginary part of the 
potential. It is thus mainly sensitive to the shape of the 
Coulomb barrier which in turn is defined by the exterior 
of the real part of the nuclear potential whereas details 
of the imaginary part are not so important. This is sim- 
ilar to the full absorption for small L at higher energies, 
but only very few a particles are able to tunnel through 
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the Coulomb barrier at very low energies. An exterior 
contribution is obtained from the tail of the imaginary 
potential for large radii; it docs practically not depend 
on the real part. This tail is weak but affects all incoming 
a particles. The relative importance of the interior and 
exterior contributions can be estimated by a variation of 
the strengths of the real and imaginary nuclear potentials 
or by introducing a cutoff radius i? C ut for the imaginary 
potential (W(r) = for r > R C ut)- 

It is found that at energies around 12 MeV an enhance- 
ment of the real part by a factor of two increases er roac 
by a factor of about 3 whereas an enhancement of the 
same factor of two for the imaginary part increases cr roac 
only by about 40%; i.e., the interior absorption is dom- 
inating at these energies. This is confirmed by the fact 
that a cutoff for the imaginary potential at i? cu t = 12 fm 
reduces cr roa c by less than 10%. 

At much lower energies the relative contribution of the 
exterior absorption increases, and at about 6 MeV it is 
found that a factor of two enhancement of the real part 
and the imaginary part both lead to the same factor of 
two enhancement for cr rea c- A cutoff radius of i? C ut = 
12fm for the imaginary part now reduces er re ac by about 
a factor of 6, and even a huge i? C ut = 15fm still leads to 
a reduction of cr rcac by almost a factor of two. 

In other words this means that the prediction of cr rea ,c 
at very low energies requires detailed knowledge of the 
nuclear potential at very large distances. The real part 
at large distances is essential for the internal contribu- 
tion at energies of about 10 — 15 MeV, and the imagi- 
nary part at large distances is important for the external 
contribution which dominates at very low energies far 
below 10 MeV. However, the sensitivity of elastic scat- 
tering data to the potential at large radii is limited. It 
is obvious that the reproduction of cr TegbC is a necessary 
requirement for any global potential; but this require- 
ment is almost automatically fulfilled at energies above 
the Coulomb barrier. Further investigation is needed on 
the question whether such a global potential is able to 
predict reaction cross sections of a-induced reactions be- 
low the Coulomb barrier. It remains an open question 
which deviation (e.g. in x 2 / F) remains allowed in the 
description of elastic scattering data to ensure a reason- 
able prediction of low-energy a-induced reaction cross 
sections. In this context an interesting example is the 
potential of Frohlich and Rauscher 3l|; it has been op- 
timized for the calculation of reaction cross sections and 
works very well for this purpose, but it shows larger de- 
viations than other global potentials in the analysis of 
elastic scattering angular distributions. 

Table IIIII shows that the relative uncertainty of cr reac 
increases strongly with decreasing energy. This is not 
surprising because the elastic scattering angular distri- 
bution approaches the Rutherford cross section at low 
energies. This holds in particular for the lowest energy 
where the ratio to the Rutherford cross section is above 
0.8 even at the most backward angles; thus, very pre- 
cise data are required here. Unfortunately, these data of 



Watson et al. [12J are not available numerically and have 
been extracted from Fig. 2 of Hi in the EXFOR Q3] 
database. Using these EXFOR data leads to poor fits 
with parameters outside the expected range and <7 rcac 
rj 65 mb. I have carefully repeated the digitization of the 
data in Fig. 2 of [lH, and I have obtained significantly 
larger elastic cross sections at backward angles. Now the 
fit provides parameters in the expected range, and I find 
a smaller total reaction cross section of a lcac = 28.1 mb. 
It is interesting to note that the original analysis in [l2[ 
using a Woods-Saxon potential gives an almost identical 
total reaction cross section of a Teac = 27.8 mb. This con- 
firms the newly digitized data. These revised data will 
be sent to EXFOR. 

For comparison of various targets at different ener- 
gies, often reduced total reaction cross sections er ro d = 

^reac 

/(^p /3 + A^ 3 ) 2 are plotted versus the reduced en- 
ergy E ied = (4/ 3 +4/ 3 )£ c . n ,/(ZpZ T ), see e.g. [HH. 
The data for many a- nucleus systems are shown in Fig. [5] 
It is obvious that the new data for 140 Co fit nicely into 
the global systcmatics. Furthermore, it can be seen that 
various global potentials [lfl [ID, [l|| provide almost iden- 
tical <7 rc d in the energy range above E le< i > 0.8 MeV 
(E c .m. 13.7 MeV) whereas significant deviations be- 
tween the different potentials are only found at very low 
energies. It should be kept in mind that at these very 
low energies (E Ted <C 0.8 MeV) the total reaction cross 
section is dominated by inelastic scattering, mainly by 
Coulomb excitation, but not by compound formation. 



C. Lower energy limit for the determination of cr re ac 

The extraction of er r cac from clastic scattering angu- 
lar distributions is possible with reasonable uncertainties 
as long as the deviation from the Rutherford cross sec- 
tion exceeds the experimental uncertainties significantly. 
The Watson data at 15 MeV fulfill this requirement with 
<j/<jr ~ 0.8 (i.e., a deviation from Rutherford of about 
20 %) at backward angles and claimed experimental un- 
certainties of less than 5%. However, at slightly lower 
energies the extraction of tr reac becomes impossible even 
with typical uncertainties of up-to-date high-precision 
scattering data [H, [3(| which are still of the order of 
a few per cent (systematic plus statistical uncertainties). 

A test calculation has been made at E\ a \, = 14.0 MeV 
(E c m . = 13.61 MeV) in the following way. First, the an- 
gular distribution at this energy is calculated from the 
ATOMKI-V1 potential from 20° to 175° in steps of 5°. 
Here I find ct/ctr ~ 0.94 at the most backward angles. 
Next, several virtual experimental data sets are created 
from this angular distribution by randomly varying each 
data point within its uncertainty (Gaussian distribution 
with an assumed 3 % Tc uncertainty, corresponding to 
the most precise available scattering data [IH, Ho|). Fi- 
nally, the various virtual data sets are analyzed in the 
usual way by fitting the parameters of the optical poten- 
tial. Here I use two parametrizations: (i) a folding poten- 
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FIG. 5: (Color online) Reduced reaction cross sections <7 rc d 
versus reduced energy £rcd for various Q-nucleus systems 
[TcL [13, [27 || . comp ared to predictions from global a-nucleus 
potentials [lfj, [isi. 19"]. The experimental data are taken from 
H2i0]l f° r 140 Ce and from [l5|] for natural cerium. The pre- 
dictions are almost identical in the shown energy range. The 
indirect data point at 10 MeV from the 143 Nd(n,a) 140 Ce reac- 
tion is located far below the shown region (E IC & = 0.584 MeV, 
ffred = 4.0 x 10~ 5 mb). 



tial in the real part and a surface Woods-Saxon potential 
in the imaginary part (i.e., the same parametrization as 
the underlying ATOMKI-V1 potential), and (ii) a vol- 
ume Woods-Saxon potential in the real and in the imagi- 
nary part. The results are shown in Fig. [6] As expected, 
it is found that the fits become relatively unstable, and 
restrictions on the number of fitted parameters have to 
be used. The obtained total reaction cross sections from 
the fits to the virtual data sets vary between about 5 and 
20 mb (i.e. vary by a factor of four) whereas the starting 
value from the ATOMKI-V1 potential was 8.8 mb. Thus, 
the uncertainty for the derived <7 rcac is at least a factor 
of two at the low energy of -E c .m. = 13.61 MeV which is 
dramatically larger than the about 15% uncertainty ob- 
tained at the slightly higher energy of 14.28 MeV. How- 
ever, the largest deviations for cr roac are correlated with 
extreme fitting parameters and unusually oscillating an- 
gular distributions (e.g. in the case of the 4 th virtual data 
set). Constraining the resulting parameters to a reason- 
able range may allow to determine cr reac even in the case 
of the 4 th virtual data set; but such constraints put into 
question whether this is still a model-independent deter- 
mination of er r0 ac • 

A simple constraint might be to fix the geometry of the 
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FIG. 6: (Color online) Rutherford normalized elastic scatter- 
ing cross sections of 140 Ce(a,a) 140 Ce versus i9 c .m. at i5 c .m. = 
13.61 MeV. Various virtual data sets are created from the 
ATOMKI-Vl potential (full black line) by a random varia- 
tion of the data points within an assumed 3 % uncertainty, 
and they are analyzed to estimate the uncertainty of the de- 
rived total reaction cross section <r rcac using either a folding 
potential in the real part (case i, blue dashed) or volume 
Woods-Saxon potentials (case ii, red dotted). Further discus- 
sion see text. 



potentials from data at higher energies, e.g. by taking the 
average values of w in the real part and Rs and as in the 
imaginary part from Table [TTJ Then only two parameters 
(A and Ws) remain to be adjusted to the data sets. How- 
ever, even in this very restricted parameter space total 
reaction cross sections a roac between about 5 and 20 mb 
have been found, and it remains necessary to keep A and 
Ws close to their expected values to obtain values of cr roac 
close to the initial cr rcac = 8.8 mb. So the lower limit for 
the determination of (T reac seems to be reached around 
E w 13.5 MeV even if high-quality scattering data are 
available over the full angular range. 
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IV. 



A SHORT NOTE ON THE 143 ND(n,a) 140 CE 
REACTION 



Recently, experimental data for the 143 Nd(n,a) 140 Ce 
reaction have been measured at energies below 1 McV 
[3^ | and between 4 and 6 MeV [33L and earlier data be- 
low lOMeV are available from [34j and [35j]. The cal- 
culation of the 143 Nd(n,a) 140 Ce reaction cross section 
in the statistical model requires the knowledge of the 
140 Ce-a potential. Because of the huge positive Q-value 
of 9.72 MeV, this potential has to be determined at en- 
ergies around 15 MeV for the higher-energy data in [33{. 
This is exactly the energy range of the Watson et al. scat- 
tering data [l2| where the various potentials under study 
were able to reproduce a reac within about 15 — 20% (see 
Table ITTT)) . Thus, the uncertainty from the a-nucleus po- 
tential in the calculation of the 143 Nd(n,a) 140 Ce reaction 
remains very limited. I have calculated the cross sec- 
tions using the default parameters of Talys VI. 4 (3(| 
and the MCF potential in the a channel. The result is 
shown in Fig. [7] The agreement with the experimental 
data above 1 MeV [HI, is reasonable although the en- 
ergy dependence of the data points is not perfectly repro- 
duced. At energies below 1 MeV the experimental data 
of [iU, HH are overestimated which is a typical behav- 
ior for the MCF potential at energies significantly be- 
low the Coulomb barrier, see e.g. for 144 Sm(a,7) 148 Gd 
Q, 141 Pr(a,n) 144 Pm d, [l|, 112 Sn(a, 7 ) 116 Te Q, and 
106 Cd(a, 7 ) 110 Cd 1. 
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FIG. 7: (Color online) Cross section of the 143 Nd(n,o^ 40 Ce 
reaction. The experimental data are taken from [32 - [35l |. 
Calculations have been performed using the MCF potential 
in the statistical model codes Talys [3y| (blue dotted) and 
Smaragd (full brown). The results for the ATOMKI- 
VI (green dashed) and AVR (red dash-dotted) potentials are 
derived from the Talys results. Further discussion see text. 



The calculations have been repeated using the code 
Smaragd [37[ with its default parameters (here the MCF 
potential is used in the a channel by default), see Fig. 
brown line. As pointed out in [H|, the 143 Nd(n,a) 140 Ce 
reaction cross section is mainly sensitive to the optical 



potential in the a channel whereas the other ingredients 
of the statistical model codes have only minor influence 
on the (n,a) cross section (for details see the sensitivity 
figures shown in Ref. [HI). As expected, the two codes 
provide very similar results for energies below 1 MeV. 
However, above 1 MeV surprisingly significant differences 
up to a factor of two appear. Interestingly, the Smaragd 
calculation shows a similar energy dependence as the ex- 
perimental data points [HI H3] but it underestimates the 
abolute scale significantly. It is beyond the scope of the 
present paper to provide a detailed comparison of the 
Talys and Smaragd codes. The disagreement between 
both codes may be related to the treatment of excited 
states in the exit channel (i.e., in the 140 Ce nucleus); 
this sensitivity is not explicitly studied in [39j | . But dif- 
ferent choices of this treatment in Talys lead only to 
minor variations of about 20 % for the (n,a) cross sec- 
tion and cannot explain the factor of two discrepancy 
between Talys and Smaragd. Therefore, great care is 
required in the determination of the a-nucleus potential 
from low-energy reaction data because the influence of 
other ingredients on the calculated cross sections has to 
be investigated in detail. 

Nevertheless, the influence of different potentials on 
the 143 Nd(n,a) 140 Ce reaction cross section can be studied 
at low energies below 1 MeV where Talys and Smaragd 
show the expected agreement for the MCF potential. It 
is found that the MCF potential overestimates the ex- 
perimental data by a factor of 5. As this energy range 
below 1 MeV corresponds to energies in the Q-channel 
of about 10 MeV, it is possible to estimate cr reac at this 
energy because cr 10 ac is only sensitive to the chosen a- 
nucleus potential: cr roac ~ 1.9 ± 0.4 [ih. For comparison, 
the predictions from the different potentials under study 
at 10 McV arc also listed in Table [TTT1 

It is technically difficult to use the ATOMKI-V1 fold- 
ing potential in the Talys or Smaragd codes. However, 
because the (n,a) cross section is practically sensitive 
only to the a-potcntial 38, 39], it is possible to calculate 
the 143 Nd(n,a) 140 Ce cross section for the ATOMKI-V1 
potential by scaling the MCF result with the ratio of the 
total cross sections a lcac of the ATOMKI-V1 potential 
and the MCF potential at the corresponding energy in 
the a channel. The same method has been applied for 
the AVR potential. In both cases the underlying MCF 
result has been taken from the Talys calculation. The 
results are compared to the experimental data of [32T4351 ] 
in Fig. [71 At energies below 1 MeV the ATOMKI-V1 
potential slightly overestimates the experimental data of 
[32l [35( 1 whereas the AVR potential slightly underesti- 
mates the data. But both potentials are much closer to 
the experimental data than the MCF potential which pre- 
dicts the (n,a) cross section about a factor of two higher 
than the ATOMKI-V1 potential and more than one order 
of magnitude higher then the AVR potential. 

For the data in the MeV region [HI, HH (correspond- 
ing to the 15 MeV region in the a channel) all potentials 
under study provide almost identical total reaction cross 
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sections cr roac , and thus the calculated (n,a) cross sec- 
tions are very similar. It is not possible to extract infor- 
mation on the potential in this energy range. This holds 
in particular as long as significant differences appear be- 
tween the two widely used codes Talys and Smaragd. 
Finally, it is interesting to note that almost the same 
behavior of the MCF (overestimation), AVR (slight un- 
derestimationVand ATOMKI-V1 (slight overestimation; 
not shown in [l(|) potentials is found in the analysis [l(| 
of the recent 141 Pr(a,n) 144 Pm data @. 

V. SUMMARY AND CONCLUSIONS 

Angular distributions of 140 Ce(a,a) 140 Ce elastic scat- 
tering were analyzed in the framework of the optical 
model. Excellent fits were obtained using locally adjusted 
parameters of the real folding and imaginary Woods- 
Saxon potentials. Total reaction cross sections er rcac were 
derived from the local fits using Eq. ([1}. The new data 
are in excellent agreement with the systematics of re- 
duced cross sections er re d [H, H3] • 

It was shown that various global a-nucleus potentials 
predict the total reaction cross sections c reac very well 
although the predicted angular distributions do not agree 
perfectly with the experimental angular distributions and 
may even deviate by up to an order of magnitude at 
very backward angles. Thus, at least at energies above 
about 15 MeV, the a-nucleus potential is not the major 
source of uncertainties in the calculation of a-induccd 
reaction cross sections in the statistical model. This can 
be understood by a semi-classical interpretation of the 
partial wave analysis. 

It is highlighted that data with very small uncertain- 
ties are required for the extraction of <7 reac at low en- 
ergies where the elastic scattering angular distribution 
approaches the Rutherford cross section. As the analy- 
sis of the Watson et al. data [12] shows, special diligence 



is required for data which are taken from published fig- 
ures by digitization. An extraction of <7 re ac from elastic 
scattering angular distributions becomes impossible with 
a reasonable uncertainty as soon as the deviation from 
the Rutherford cross section at backward angles becomes 
smaller than the experimental uncertainty. 

Finally, further information on the a-nucleus poten- 
tials at very low energies can be extracted only from 
the analysis of reaction data. Although great care is re- 
quired for such an analysis because of the influence of 
other ingredients of the statistical model calculations, it 
can be concluded from the low-energy 143 Nd(n,a) 140 Ce 
data below E n = 1 MeV that the recent ATOMKI-V1 
and AVR potentials provide a significant improvement 
around E ps 10 MeV in the a channel, i.e. far below the 
Coulomb barrier, whereas the MCF potential overesti- 
mates the 143 Nd(n,a) 140 Ce data by a factor of 5. 
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